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ABSTRACT 

The observed radial velocity (RV) eccentricity distribution for extrasolar plan- 
ets in single-planet systems shows that a significant fraction of planets are ec- 
centric (e > 0.1). However an RV planet's eccentricity, which comes from the 
Keplerian fitting, can be biased by low signal-to-noise and poor sampling. Here 
we investigate the effects on eccentricity produced by undetected outer compan- 
ions. We have carried out Monte Carlo simulations of mock RV data to under- 
stand this effect and predict its impact on the observed distribution. We first 
quantify the statistical bias of known RV planets' eccentricities produced by un- 
detected zero-eccentricity wide-separation companions and show that this effect 
alone cannot explain the observed distribution. We then modify the simulations 
to consist of two populations, one of zero-eccentricity planets in double-planet 
systems and the other of single planets drawn from an eccentric distribution. Our 
simulations show that a good fit to the observed distribution is obtained with 
45% zero-eccentricity double-planets and 55% single eccentric planets. Assuming 
our two simulated populations of planets are a good approximation for the true 
RV population, matching the observed distribution allows us to determine the 
probability that a known RV planet's orbital eccentricity has been biased by an 
undetected wide-separation companion. Averaged over eccentricity we calculate 
this probability to be ~ 4%, suggesting that a small fraction of systems may have 
a yet to be discovered outer companion. Our simulations show that moderately- 
eccentric planets, with 0.1 < e < 0.3 and 0.1 < e < 0.2, have a ~ 13% and 
~ 19% probability, respectively, of having an undetected outer companion. We 
encourage both high-contrast direct imaging and RV follow-up surveys of known 
RV planets with moderate eccentricities to test our predictions and look for pre- 
viously undetected outer companions. 



Subject headings: planetary systems: radial velocities-eccentricity 
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1. Introduction 

Since the 1990s, over 300 extrasolar planets have been discovered (http:/ /exoplanet.eu). 
Exoplanets have been discovered by several different techniques, namely radial velocity (RV), 
transits, astrometry, microlensing, and most recently direct imaging (e.g., Kalas et al. (2008); 
Marois et al. (2008)). The RV method, which has discovered the majority, monitors the 
periodic velocity shifts (in the radial direction) of a given star. These velocity variations 
are caused by the star orbiting the barycenter of the combined system, which can contain 
additional companions (typically stars, brown dwarfs, and/or planets). The RV method 
allows one to extract the orbital period of the companion (P), its semi-major axis a, the 
lower limit to its mass (mshu), and its eccentricity (e). Since there are ~ 300 RV-discovered 
exoplanets, we have a large enough sample to construct distribution functions for each of 
these quantities. 

Of these distributions, the eccentricity distribution is the most puzzling. In our own 
solar system, planets have very low eccentricity values (nearly circular orbits). Jupiter, the 
gas giant to which we compare discovered extrasolar planets, orbits with e = 0.05. While 
the RV extrasolar planet eccentricity distribution has a strong peak near e = 0, it also has 
a significant tail extending all the way out to e = 0.93 (see Fig. 1). If our solar system and 
its planets are common and ordinary, then we would expect most exoplanets to have low 
eccentricities. Clearly this is not the case. 

There have been many attempts to explain the observed eccentricity distribution. Most 
of these focus on planet-planet interactions and scattering (e.g. Rasio & Ford (1996); Levison 
et al. (1998); Adams & Laughlin (2003); Juric & Tremaine (2008)). These interactions can 
build up the eccentricities of two or more planets, resulting in scattering or ejection. In 
the process, the remaining planets settle down into stable, high-eccentricity orbits. While 
scattering is certainly plausible and can match the e > 0.2 observed eccentricity distribution, 
it is unclear whether the majority of planetary systems go through such violent planet-planet 
interactions. 

Other factors may be at play. Shen & Turner (2008) explored the effects of RV sys- 
tematics on eccentricity. They showed that planets' eccentricities can increase with lower 
signal-to-noise ratios (SNR), given by K/a, where K is the velocity semi-amplitude and a 
is the experimental error (typically depending on the telescope and instrument error as well 
as stellar jitter). Specifically e increases for K/a < 3. They also show that eccentricity 
increases with fewer observations in a given RV data set (N obs < 60). Here, the number of 
observations is the number of data points in a star's radial velocity curve. Since the SNR 
and N obs conditions are not always met for RV data, especially for the earliest discoveries, 
we can expect that this bias plays a role in the observed eccentricity distribution. Shen & 
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Turner (2008) estimated that about 10% of RV-detected exoplanets are affected by this bias. 

Undetected additional planets in planetary systems can also affect RV eccentricities. 
When two planets orbit a star, both contribute to the star's radial velocity. Assuming the 
planets are not gravitationally interacting with each other, one produces the total radial 
velocity curve by adding the radial velocity data from each planet. One can detect the 
second planet by monitoring the star long enough to observe its full period. Since this time 
frame is usually > 10 years, only ~ 20 multiple-planet systems have been discovered. If one 
observes a double-planet system for a duration shorter than the period of the second planet, 
a long-term trend can appear in the data. Fischer et al. (2001) explored the possible bias 
introduced by wide- separation planets by injecting real long-term trend data into some of 
their existing one-planet RV data. After fitting the new data, they found that the mass and 
semi-major axis of the original planet were mostly unchanged, but that its eccentricity value 
increased, sometimes by as much as Ae = 0.25. They were able to detect the long-term trends 
only in cases where there was dense phase coverage (large N f, s ) and high SNR. These results 
are very important for several reasons: they confirm the effects of the systematic biases 
discussed by Shen & Turner (2008), and they show that undetected additional companions 
can drive up a known planet's eccentricity value. Nonetheless, we cannot draw large-scale 
conclusions since the Fischer et al. (2001) results were based on a small sample of tests, and 
the test parameters were not significantly (randomly) varied. 

To determine the full effect of undetected wide-separation planets on eccentricity, we 
have carried out Monte Carlo simulations of mock RV data. We first determine the statistical 
effect of undetected wide-separation companions on eccentricity. We then try to match the 
observed eccentricity distribution with a population of zero-eccentricity planets in double- 
planet systems. We show that this population alone cannot describe the observed eccentricity 
distribution, and we therefore need a second input population of high-eccentricity planets. 
After matching the observed distribution, we derive the probability that a known planet has a 
wide-separation companion, information potentially important for both RV follow-up surveys 
and direct imaging. In §2, we describe the simulations, test parameters, and methodology. 
In §3, we present the results and discuss the implications. In §4, we summarize and conclude. 



An RV-detected planet's orbital parameters (msini, a, and e) are determined from the 
Keplerian fit to the star's radial velocity data. The equations needed to calculate the star's 
radial velocity are: 



2. Simulation Parameters and Methodology 




(1) 
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K = yjG/a/{M* + m) m sin i / y/l - e 2 



(2) 
(3) 
(4) 



T = 2 tan _1 [tan(E/2) + e)/(l - e)] 



£ = M + e sin £ 



M =(*-«„) — . 



(5) 



P is the orbital period of the planet derived from Kepler's laws, t is the time of the obser- 
vation, t is the time of perihelion, and M is the mean anomaly. E is the eccentric anomaly 
and since it is given in a transcendental equation it cannot be solved analytically. For our 
simulations, we solve for E using iterative loops. T is the true anomaly, K is the velocity 
semi-amplitude, and uo is the argument of the perihelion. We have set sin? = 1, assuming 
edge-on systems in all cases. Vq is the velocity offset which for simplicity we have set to 
zero. M* and m are the mass of the star and planet, respectively. Since most RV stars are 
solar- type stars, we set M* = M Q . 

Our simulations, aside from control simulations, consist of 100 sets of N RV curves. N = 
227 since, as of January 2009, there are 227 single-planet RV systems (http:/ /exoplanet.eu). 
Thus we simulate ~ 10 4 total planets. We do not simulate the total number of observed RV 
systems because some of these (~ 30 systems) have multiple planets. Our hypothesis applies 
only to single-planet RV detections. 

For each simulated planet, we draw oo randomly from a uniform distribution in the range 
[0,27r] and to i n the range [0,P], where P is the period of that particular planet. Each planet's 
semi-major axis a and mass m are drawn randomly from the most current RV distribution 
functions, irregardless of being in a single or double-planet system. From the Butler et al. 
(2006) catalog of exoplanets, these are dN/dM oc M _L1 and dN/d\oga oc a 0A . In simulated 
systems, the first planet has ai nner G [0.05, 3] AU and the second has a outer G [7, 15] AU. 
The first planet's orbital range was chosen because most RV-detected planets are found 
within 3 AU, and the closest planets have periods on the order of a few days (= 0.05 AU). 
We acknowledge that the Butler et al. (2006) minimum semi-major axis used to derive the 
distribution function was a ~ 0.1, but this difference is insignificant. It is also true that 
RV has found some planets beyond 3 AU, but the statistics are poor for this small sample. 
The second planet's orbital range was chosen for several reasons: gas giants in our solar 
system (Jupiter, Saturn) orbit in this range, and we wanted to avoid any resonances < 2:1 
that might, via secular interactions, lead to higher eccentricity values. In fact, the smallest 
resonance we allow is ~ 3.5 : 1 (outer planet at 7 AU, inner planet at 3 AU). Our simulations 
are testing stable relaxed planetary systems, and the eccentricity bias we are investigating 
is not from planet-planet interactions but rather from the Keplerian fitting. 
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The mass range for the first planet was chosen to be m inner E [0.01, 25] Mj (Mj is the 
mass of Jupiter). These limits reflect the ~ minimum and maximum mass planets discovered 
by RV so far (http://exoplanet.eu). The mass range for the second planet was chosen to be 
m outer £ [ m inner/Mj, 25] Mj. If a planetary system has two planets, then the outer planet 
is probably at least as massive as the inner. Though speculative, there is some evidence for 
this in Wright et al. (2008), who show that there is on average a 1-1 correspondence in mass 
for the RV-discovered two-planet systems. It is possible that this result is biased by the 
fact that more distant planets have to be more massive to be detected by the RV method. 
Nonetheless most multiple-planet systems have been discovered by RV, so the sample in 
Wright et al. (2008) is the largest and most statistically-significant sample from which to 
draw conclusions. 

For each simulation, we compute the star's radial velocity V ra d at 32 points in time 
over 2P. In all cases (one-planet system or two-planet system), P is the period of the inner 
planet. We have chosen 32 data points because this is a moderate amount of observations. We 
acknowledge that since N Q b s = 32 < 60, eccentricities will be slightly biased (as described 
by Shen & Turner (2008)). Nonetheless, N Q b s = 32 is a reasonable assumption given the 
typical N b s for an actual RV data set. We have chosen to simulate each system's data for 
two periods because most RV data is published with at least two periods of coverage. We 
emphasize that the 32 data points along each RV curve are not equally spaced. This would 
not be realistic. Instead we have tried to mimic the "clustering" of data points seen in 
typical RV data (see Fig. 2). To accomplish this, we divide a given simulated data set into 
8 zones, each containing 4 points (observations). In each zone, the first point's location is 
determined randomly. The subsequent points' locations are also determined randomly, but 
the inner boundary for each is the location of the previous point. This causes an apparent 
grouping, or clustering, of observations in each zone. 

If there are two planets in a simulated system, we compute the radial velocity inde- 
pendently for each planet and then add them together to get the total V ra d- In all cases, 
we add Gaussian noise with a — 5 m/s to V ra d- This a was chosen as a good estimate of 
typical RV precision. Because we enforce K > 3cr, changing a to a higher (potentially more 
RV-representative) value (say, 10 m/s) would have little to no effect on our results. All sim- 
ulations were carried out in Matlab. To determine the orbital parameters, we cycle through 
e from to 0.99, hold e constant at each value, and fit for the other orbital parameters. We 
then take the lowest x 2 solution, along with its corresponding e, as the final set of orbital 
parameters. Because we know the input parameters for each simulated system beforehand, 
we do not use periodograms and false-alarm probability tests in our fitting procedure. When 
fitting a data set, we only fit for one planet because we are testing the effects of undetected 
wide-separation planets on eccentricities. To ensure that the outer planets do not create 



- 6- 



long-term trends in the data, we enforce Xreduced < 3- This limiting value was chosen be- 
cause RV data with poorly-constrained fits are rarely published; if they are, the authors 
typically suggest that a long-term trend is present and attempt a two-planet Keplerian fit. 

3. Results and Discussion 

3.1. Control Simulations 

Before running our full suite of simulations, we checked that our Keplerian fitting was 
working correctly. To test this, we simulated 50 sets of N = 227 single-planet systems (~ 
10,000 planets). Each planet's semi-major axis and mass were drawn from the appropriate 
distributions. We set e to zero, which allows us to directly measure any unforeseen effects 
on eccentricity. If the fit is working correctly, then the output parameters for each planet 
should match the input parameters within experimental error. 

Fig. 3 shows the relationships between a ou t pu t and ai nput , and m output and m inpu t. The 
tight linear relationships and slopes of unity show that our Keplerian fitting has no effect or 
bias on a planet's semi-major axis or mass. Fig. 4 shows the average distribution of fitted 
eccentricity values, normalized to the first bin. Fig. 4 reveals that ~ 8% of planets are 
eccentric (e > 0.1). This small bias is due to the effects discussed by Shen & Turner (2008). 
Namely, our N b s = 32 < 60, which is the suggested number of observations. Furthermore, 
although Shen & Turner (2008) did not explore this, we believe that the "clustering" of data 
points also affects a planet's eccentricity. This is a typical feature of published RV data, so 
we have chosen not to ignore it. 

Having confirmed the validity of our Keplerian fitting routine, the next step was to 
quantify the statistical effects of undetected wide-separation planets on one-planet fits. To 
this end, we simulated 100 sets of N ~ 100 two-planet systems (~ 10, 000 planetary systems). 
Each planet's semi-major axis and mass were drawn from the appropriate distributions, 
and we again fixed e at zero. Comparing the results of these double-planet simulations to 
those of the single-planet simulations yields the second planet's effects on the first planet's 
parameters. Fig. 5 shows the relationships between a output and a input , and m output and m input . 
In all cases, a and m are the semi-major axis and mass of the inner planet. Fig. 5 shows tight 
linear relationships with slopes of unity. Evidently wide-separation planets have no effect on 
the inner planet's mass and semi-major axis, although there is certainly more scatter. These 
results agree with the previous work by Fischer et al. (2001). 

Fig. 4 also shows the average distribution of fitted eccentricity values for two-planet 
systems. The percentage of planets with 0.1 < e < 0.2 is ~ 13%, and the total percentage 
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of planets with e > 0.1 is ~ 18%. Thus the addition of the second planet can increase the 
eccentricity of the first planet. Knowing that the N b s < 60 effect accounts for ~ 8% of the 
biased eccentricities, the sole effect from wide- separation planets is ~ 10%. We note that 
"wide-separation" means 7 < a outer /A\J < 15, the final range we settled on. We have also 
tested other semi-major axis ranges of different extents, such as [10,30] AU. In general, we 
saw that the larger the semi-major axis range, the lower the eccentricity bias. This makes 
sense given the RV equations (Eqs. 1-5); closer planets increase the velocity semi-amplitude 
K which in turn affects V ra d and the parameters of the fitted planet. In our simulations, 
however, bringing the second planet closer in does not necessarily induce more eccentricity 
bias. This is because the Keplerian fit becomes poor and the system is rejected. Therefore 
there is a "sweet spot" (~ 7 AU) where eccentricity is slightly biased, no long-term trend 
appears in the data, and the Keplerian fit converges. 

The double-planet eccentricity distribution in Fig. 4 led us to draw two additional con- 
clusions. First, our results agree with those of Fischer et al. (2001). The particular scatter 
in eccentricity values (e as high as ~ 0.3) matches the maximum Ae they reported. If we 
relaxed our fitting constraints, then we would surely see higher eccentricity values. Therefore 
our agreement with the Fischer et al. (2001) result suggests that our fitting constraints (in 
particular Xreduced < 3) are acceptable. Second, we cannot match the observed single-planet 
RV eccentricity distribution (Fig. 1) with only zero-eccentricity double-planet systems. Fig. 
4 at best only matches the e < 0.2 observed distribution shown in Fig. 1. Consequently 
we can reject the notion that the observed distribution can be explained by 100% zero- 
eccentricity double-planet systems. This result has two interesting implications. It lends 
more support to the growing evidence that the true RV population consists of a significant 
fraction of high-eccentricity planets in single-planet systems; and, following from this con- 
clusion, our solar system, as current evidence continues to show, is probably not ordinary. 
This is to say that the number of planets in our solar system and the low eccentricity of 
each is uncommon. If we had been able to match the observed eccentricity distribution with 
227 zero-eccentricity double-planet systems, we could suggest that multiple-planet systems 
with low-eccentricity planets (like our own) may be more common than previously believed. 
However our simulations show that multiple-planet systems cannot be the entire story. 



3.2. Matching the Observed Eccentricity Distribution 

Our preliminary results show that we cannot match the observed single-planet eccentric- 
ity distribution with only zero-eccentricity double-planet systems. High-eccentricity planets 
in single-planet systems must constitute a significant fraction of the true RV population. 
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Therefore, in addition to a simulated population of zero-eccentricity double-planets, we intro- 
duce a second population of single truly-eccentric planets. When simulating a truly-eccentric 
single planet, we draw its eccentricity value randomly from the Schwarzschild distribution 
for eccentricity e, given by 

2 

dN/de=^e W (^), (6) 

with a e = 0.3. This is the distribution in eccentricity one would expect to see from gravi- 
tational scattering of astronomical bodies, whether they are stars or planets. Furthermore, 
Juric & Tremaine (2008) used it to generate eccentricities in their simulations, and it matches 
the e > 0.2 observed eccentricity distribution very well. 

Fig. 6 shows the results of our final simulations, comparing three simulated distributions 
to the observed single-planet eccentricity distribution from http://exoplanet.eu. The ratio 
of double-planet systems to total systems (= R) varies for each of the different simulations 
between 30%, 45%, and 65%. Fig. 6 shows that the R = 30% case underestimates the number 
of low-eccentricity planets and overestimates the number of high-eccentricity planets. The 
R = 65% case shows the opposite. The R = 45% distribution is visibly the best match, with 
each point at or within 1 sigma of the observed distribution. This case also had the lowest 
Xreduced value (= 1.74) compared to the other two distributions, which had Xreduced > 2. Fig. 
7 shows the cumulative distribution functions of the simulated and observed distributions. 
Again the overall agreement is evident for the R = 45% distribution. We also performed 
a two-sample K-S test for the R = 45% case. The null hypothesis was that this simulated 
distribution and the observed distribution from http://exoplanet.eu were the same. The 
two-sample K-S test showed that the null hypothesis could not be rejected at the a = 0.05 
significance level (95% confidence), supporting the notion that the R = 45% distribution 
matches the observed distribution very well. These statistical results, as well as the final 
simulation parameters, are detailed in Table 1. 

To zero in on the optimal ratio R, we carried out several simulations that varied R 
around 45%. Together with the R = 30% and 65% simulated distributions, we were able 
to quadratically fit Xreduced as a function of R. This yielded a minimum Xreduced °f 1-6, 
corresponding to R m %n = 42+®;^% (1-sigma limits). 

In a two-planet system, what determines if the fitted planet's eccentricity is biased? 
We already know the effect depends on the semi-major axis and mass of the second planet 
from Eqs. 1-5, but which of these is the dominating factor? Fig. 8 shows the ratio of the 
semi-major axes and the ratio of the masses of the two planets in each double system as 
a function of fitted eccentricity. From Fig. 8 it is evident that a lower e means a larger 
separation between the two planets, and a higher e corresponds to a smaller separation. 
There is no apparent correlation in mass with eccentricity. This tells us that in a one-planet 
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Keplerian fit, the fitted eccentricity value depends much more on how far away the second 
planet is than on how massive it is. 

3.3. Eccentricity as an Indicator of Multiplicity 

The ratio R tells us the fraction of simulated systems that have wide-separation com- 
panions, but we do not know which known RV systems might be affected by outer-planet 
bias. Because we have succeeded in matching the observed eccentricity distribution, we can 
use eccentricity as an indicator of exoplanet multiplicity. We stress that the discussion here- 
after assumes that our two simulated input distributions of planets are a good approximation 
for the true RV population. To determine what information about multiplicity a planet's 
eccentricity yields, we calculate the probability that an RV-detected planet's eccentricity has 
been biased by a wide-separation companion. Fig. 9 plots the ratio of the average number 
of double-planet systems to total systems as a function of fitted eccentricity. The first bin 
(0 < e < 0.1) has a possible degeneracy in eccentricity: a planet with an eccentricity in this 
range either lives alone or has an extremely distant companion. In both cases, the eccentric- 
ity value is low and unbiased. We have therefore set the first bin to zero by subtracting off 
the control simulation bias (see Fig. 4). In this way any unforeseen bias introduced in our 
simulations would not corrupt our probabilities. 

From Fig. 9, the probability that a given planet with any eccentricity (0 < e < 1) has 
an undetected outer companion is ~ 4%. For e > 0.1, any wide-separation companions have 
small enough orbital separations to induce effects on fitted eccentricity. For 0.1 < e < 0.4, 
the average probability is ~ 10%. This value increases to 13% and 19% if we consider planets 
with 0.1 < e < 0.3 and 0.1 < e < 0.2, respectively. 

Thus we now have a specific range in eccentricity values that can be used for observation 
target selection. From http://exoplanet.eu, there are 84 planets (in single-planet systems) 
with 0.1 < e < 0.4, 57 planets with 0.1 < e < 0.3, and 30 planets with 0.1 < e < 
0.2. Therefore, with the corollary stated above, we estimate that about 16 known RV 
planets should have wide-separation companions. We encourage high-contrast direct imaging 
observations of these targets to look for potential wide-separation companions. For a small 
target sample size, we suggest that direct imaging surveys focus on RV planets with 0.1 < 
e < 0.2 or 0.1 < e < 0.3 where the probabilities are highest. In this case, our simulations 
predict that between 6-10 RV planets should have wide-separation companions, depending 
on which eccentricity range is used. 

Our results are important for RV follow-up surveys as well. We again suggest looking at 
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planets with 0.1 < e < 0.3. For the zero-eccentricity planets, there is certainly the possibility 
that they are not part of multiple-planet systems, in which case RV follow-up would yield 
null-results. If these planets do have wide-separation companions, they are probably so far 
away that RV follow-up would require extremely long observation time frames to detect 
anything. However for e > 0.1, RV follow-up could reveal long-term trends and potentially 
extract the periods of wide-separation companions, since our simulations predict that such 
planets would not be too distant relative to the inner planets (see Fig. 8). We summarize 
the suggested observational strategy for RV follow-up and direct imaging in Table 2. 



3.4. Why Not Observe Zero-Eccentricity RV Planets? 

One might wonder why we do not advocate direct imaging of known zero-eccentricity 
RV planets. After all, Fig. 8 shows that the two-planet systems with the most distant 
companions correspond to inner planets with nearly zero eccentricity. The more distant a 
planet is from its host star, the easier it is to directly image. However, because the outer 
planets have no effect on the inner planets' eccentricities, we cannot obtain any information 
about the systems' multiplicities; we have no way to determine if a zero-eccentricity RV 
planet lives alone or has a distant companion. Therefore if one had to choose between 
directly imaging moderately-eccentric RV planets (where we can quantify the probability 
that it has an outer companion) or zero-eccentricity planets, we recommend observing the 
moderately-eccentric planets. We cannot rule out the possibility that the RV population 
consists of a separate population of zero-eccentricity planets in single-planet systems. We 
do not entirely discourage direct imaging of zero-eccentricity planets, but the statistically 
favorable choice is to observe moderately-eccentric planets (see Table 2). 



3.5. The Uniqueness Problem 

Our simulations and results make interesting predictions, but the distributions used 
to fit the observed data are not unique. It is possible that other models with different 
parameters and assumptions could similarly reproduce the observations. Nonetheless we 
feel our results are important and worth further investigation because we have approached 
the problem from observational, theoretical, and physical standpoints. The distribution 
functions for semi-major axis and mass were obtained from statistical analysis of observed 
RV data. The orbital ranges and masses were determined from both observational data and 
theoretical assumptions. These were: in stable relaxed planetary systems, planets should 
orbit outside 2:1 resonances; and if a planetary system has two planets, the outer should be 
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at least as massive as the inner (Wright et al. 2008). Furthermore, with time the number of 
multiple-planet systems (and systems with long-term trends) discovered by RV has continued 
to increase, hinting that multiple-planet systems may be more common than the current 
statistics show. From the physical standpoint, if not all planetary systems have significant 
planet-planet interactions and scattering, then we might reasonably expect there to be two 
independent distributions of planetary systems. The planets unaffected by scattering would 
have e ~ 0, while the scattered population of planets would have eccentricities drawn from 
the Schwarzschild distribution as we have modeled. In this case our two distributions of 
planets provides a good description of the RV planet population. The predicted frequency 
of an outer companion provides the key observational test for this conclusion. 



4. Conclusions 

We have carried out Monte Carlo simulations of mock RV data to determine the statis- 
tical effect of undetected wide-separation companions on eccentricity. We quantify this effect 
and show that, for a given simulated population of planets, the number of eccentric planets 
increases by ~ 10% when the outer companion is introduced. Thus undetected outer com- 
panions in RV planetary systems can have a small but important effect on the inner planet's 
Keplerian-fitted eccentricity. We show that the observed RV single-planet eccentricity dis- 
tribution cannot be matched by 100% zero-eccentricity planets in double-planet systems, 
but rather an additional population of high-eccentricity planets in single-planet systems is 
required. This result lends support to the growing evidence that our solar system, with 
many low-eccentricity planets, is probably not ordinary. It also agrees with the current RV 
statistics, which show that high-eccentricity planets constitute a significant portion of the RV 
population. This is supported by our simulations requiring low-eccentricity double-planets 
to constitute less than half of the total population, with the best-fitting ratio R = 42t.f%%- 

Assuming our two simulated populations of planets are a good approximation for the 
true RV population, we can make predictions about RV planet multiplicity using eccentricity. 
To do this we calculate the probability that a known RV planet's eccentricity has been biased 
by a wide-separation companion. Averaged over eccentricity this probability is ~ 4%. For 
0.1 < e < 0.3 and 0.1 < e < 0.2, the probability is 13% and 19%, respectively. To test our 
predictions about the true RV population of planets, RV exoplanet multiplicity, eccentricity 
as an indicator of multiplicity, and to look for previously undiscovered outer companions, we 
encourage both high-contrast direct imaging and RV follow-up surveys of known RV single- 
planets with moderate eccentricities. 



- 12 - 



We would like to thank Matthew Kenworthy for helpful suggestions, comments, and his 
careful proofreading of this manuscript. We would also like to thank Jared Males and Andy 
Skemer for helpful insights into radial velocity methods and statistical analysis. We thank the 
anonymous referee for very helpful comments and suggestions which greatly improved this 
paper. This work has been supported by the NASA Astrobiology Institute under Cooperative 
Agreement No. CAN-02-OSS-02 issued through the Office of Space Science. 

REFERENCES 

Adams, F. C. & Laughlin, G. 2003, Icarus, 163, 290 
Butler, R. P., Wright, J. T., Marcy, G. W., et al. 2006, ApJ, 646, 505 
Fischer, D. A., Marcy, G. W., Butler, R. P., et al. 2001, ApJ, 551, 1107 
Juric, M. & Tremaine, S. 2008, ApJ, 686, 603 

Kalas, P., Graham, J. R., Chiang, E., et al. 2008, Science, 322, 1345 
Levison, H. F., Lissauer, J. J., & Duncan, M. J. 1998, AJ, 116, 1998 
Marois, C, Macintosh, B., Barman, T., et al. 2008, ArXiv e-prints 
Rasio, F. A. & Ford, E. B. 1996, Science, 274, 954 
Shen, Y. & Turner, E. L. 2008, ApJ, 685, 553 

Wright, J. T., Upadhyay, S., Marcy, G. W., et al. 2008, ArXiv e-prints 



This preprint was prepared with the AAS IATgX macros v5.0. 



-13- 



TABLE 1 

Summary of simulation parameters and results 



Number of planetary systems (N) 227 

Number of simulations 100 

m inner [0.01,25] Mj 

TFlouter \j^inneri 25 Mj\ 
dinner 

[0.05,3] AU 

Clouter [7,15] AU 

(^single [0)1] 

^double 

Ratio of double-planet systems to total systems (R) 45% 

^(■reduced ^ ' 

Interpolated R min 42±l%% 

Interpolated minimum Xreduced 1-6 
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TABLE 2 



RV follow-up / direct imaging observing strategy 



RV planet's eccentricity 


Comments 


RV Follow-up 


Direct Imaging 




eccentricity yields no 






e < 0.1 


information on likelihood 
of outer companion 


NO 


? 




~ 10-20% chance of 






0.1 < e < 0.4 


having an outer 
companion 


YES 


YES 




extremely unlikely 






0.4 < e < 1 


to have outer 
companion 


NO 


NO 



- 15 - 



120 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

eccentricity 



Fig. 1. — The observed RV eccentricity distribution (black bars) and the observed single- 
planet RV eccentricity distribution (white bars), both from http: / /exoplanet.eu as of January 
2009. We compare our simulations only to the single-planet distribution because we hypothe- 
size that undetected outer companions bias one-planet Keplerian fits. Note the large number 
of high eccentricity (e > 0.1) planets. 
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50 rr 




Time (days) 

Fig. 2. — Typical radial velocity plot from the simulations. In this case there are two planets 
in the system, both contributing to the star's radial velocity. There is no apparent long-term 
trend from the second planet evident in the data, but the fitted eccentricity has increased 
from e = 0toe = 0.11. 
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Fig. 3. — Left: output semi-major axis as a function of the input semi-major axis for the 
control simulation (one-planet systems). Right: the same, except showing the relationship 
between output and input mass. In each case, the tight linear relationships and slopes of 
unity verify that our simulations have no unforeseen effects on a given planet's semi-major 
axis or mass. 
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Fig. 4. — Eccentricity distribution for the control simulation (one-planet systems) (white 
bars), and the eccentricity distribution for simulated two-planet systems (black bars), both 
normalized to the respective first bin. All planets had their input eccentricities set to zero. 
For one-planet systems, 8% of the simulated planets' fitted eccentricities have increased, most 
likely due to biases coming from low SNR and poor sampling (discussed in Shen & Turner 
(2008)). When the wide-separation companions are introduced, the fraction of eccentric 
(e > 0.1) planets increases by ~ 10% to ~ 18%. 
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Fig. 5. — Left: output semi-major axis as a function of the input semi-major axis in simulated 
systems with two planets, both with zero eccentricity. Right: the same, except showing the 
relationship between output and input mass. In both cases, some outliers have been clipped. 
The tight linear relationships and slopes of unity show that wide-separation companions have 
a negligible effect on the fitted semi-major axis and mass, which agrees with expectations and 
the previous work by Fischer et al. (2001). It is worth noting, however, the increased scatter 
in output semi-major axis as compared to output mass. This suggests that semi-major axis 
is the more easily affected quantity. 



-20- 



150 



150 



* 



CD 

| 100 

Z5 



50 





A 


ft 


* 




* 












i 



-T%~U r 



0.5 

eccentricity 



CD 
.Q 

E 

Z5 



CD 
D) 
CO 

1 

CD 
> 
< 



100 



50 



0.5 

eccentricity 



150 




Fig. 6. — Simulated eccentricity distributions (skinny white bars with markers and error bars 
in all three plots) and the observed single-planet eccentricity distribution (wide grey bars 
in all three plots) from http://exoplanet.eu. Top-left: simulated distribution for R = 30%. 
The number of low-eccentricity planets is underestimated and the number of high-eccentricity 
planets is overestimated. Top-right: simulated distribution for R = 65%. This time the num- 
ber of low-eccentricity planets is overestimated and the number of high-eccentricity planets is 
underestimated. Bottom: simulated distribution for R = 45%. This distribution matches the 
observed single-planet eccentricity distribution very well (Xreduced = 1-74, disagreeing only 
around e ~ 0.55 where the observed distribution dips down, most-likely due to small-number 
statistics). 
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Fig. 7. — Cumulative distribution functions of the simulated eccentricity distributions for 
R = 30% (dotted line), 45% (dashed line), 65% (dot-dashed line), and the observed single- 
planet eccentricity distribution (solid line). The R = 45% curve matches the observed, but 
the R = 30% and R = 65% do not agree. 
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Fig. 8. — Left: ratio of the semi-major axes of the two planets in each simulated two-planet 
system as a function of fitted eccentricity. Right: the same, except the ratio of the masses. 
The semi-major axis ratio increases for low eccentricity and decreases for high eccentricity, 
but there is no apparent correlation between mass ratio and eccentricity. This suggests that 
fitted eccentricity is much more dependent on the separation between two planets than on 
the ratio of their masses. 
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Fig. 9. — Left: ratio of double-planet systems to total systems as a function of fitted eccen- 
tricity. We cannot take these results at face value because the RV systematic biases (low 
SNR, poor sampling) have not been accounted for and subtracted off. Right: probability 
that a known RV planet's eccentricity has been biased by an undetected wide-separation 
companion. To produce this plot, the control simulation bias (see Fig. 4) has been nor- 
malized to the first bin and subtracted off, explaining the first bin's zero probability. In 
actuality planets with < e < 0.1 have some small non-zero probability, but the degeneracy 
in eccentricity (single-planet system with e ~ or double-planet system with e ~ and 
an extremely distant companion) prevents us from drawing any conclusions about double- 
planet probability. For e > 0.1, there is a non-zero probability that a given RV planet has a 
wide-separation companion. Specifically, for 0.1 < e < 0.4 the average probability is ~ 10%, 
for 0.1 < e < 0.3 the average probability is ~ 13%, and for 0.1 < e < 0.2 the probability is 
~ 19%. Averaged over all eccentricities, the probability is ~ 4%. 



